ALGEBRO-GEOMETRIC SOLUTION 
OF THE DISCRETE KP EQUATION OVER A FINITE FIELD 
OUT OF A HYPERELLIPTIC CURVE 



MARIUSZ BIALECKI AND ADAM DOLIWA 

Abstract. We transfer the algebro-geometric method of construction of solutions 
of the discrete KP equation to the finite field case. We emphasize role of the Jacobian 
of the underlying algebraic curve in construction of the solutions. We illustrate in 
detail the procedure on example of a hyperelliptic curve. 



1. Introduction 

Cellular automata are dynamical systems on a lattice with values being discrete 
(usually finite) as well. They are one of the more popular and distinctive classes of 
models of complex systems. Introduced in various contexts [2B1I2S1 around 1950 they 
have found wide applications in different areas, from physics through chemistry and 
biology to social sciences [29] . 

One of the most interesting properties of cellular automata is that complex patterns 
can emerge from very simple uptade rules. However, usually one cannot easily predict 
how a given cellular automaton will behave without going through a number of time 
steps on a computer. Due to their completely discrete nature, cellular automata are 
naturally suitable for computer simulations, but also here it would be instructive to 
have examples of rules with large classes of analytical solutions, integrals of motions 
and other "integrable features". 

The problem of construction of integrable cellular automata is not new and was 
undertaken in a number of papers (see, for example, El 13 I.26J). In particular 
in [23 [TBI it was given a systematic method, called ultra-discretization, of obtainig 
cellular automaton version of a given discrete integrable system. 

Recently a new approach to integrable cellular automata was proposed in 0. Its 
main idea is to keep the form of a given integrable discrete system but to transfer the 
algebro-geometric method of construction of its solutions [THH] from the complex field 
C to a finite field. This method, which in principle can be applied to any integrable 
discrete system with known algebro-geometric method of solution, has been applied to 
the fully discrete 2D Toda system (the Hirota equation) and in [H] to discrete KP and 
KdV equations (in Hirota form). In particular, the finite field valued multisoliton 
solutions of these equations have been constructed. We remark that the algebraic 
geometry over finite fields, although conceptually similar to that over the field of 
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complex numbers [10], has its own tools and peculiarities [El EI]- It is also nowadays 
very important in practical use in modern approaches to public key cryptography [TT?] 
and in the theory of error correcting codes |25| . 

The aim of this paper is to study in the finite field context the very distinguished 
example of integrable system — the discrete KP equation. We present the algebro- 
geometric scheme of construction of its solutions in a finite field and we demonstrate 
its linearization on the level of the abstract Jacobi variety of the corresponding al- 
gebraic curve. We illustrate details of the construction on example of a hyperelliptic 
curve. 

We remark that in [3| it was observed that the Lax representation of the discrete 
sine-Gordon equation of Hirota [H] has a meaning also when the field of complex 
numbers is replaced by a finite field of the form F p 2 , where p is a prime number. The 
authors of [3] showed also that, in principle, the corresponding integrals of motion 
can be calculated. Finally, we note that possibility of considering the soliton theory 
in positive characteristic has been anticipated in |22j . 

The organization of the paper is as follows. In Section[2]we first summarize the finite 
field version of Krichever's construction of solutions of the discrete KP equation, then 
we present its abstract Jacobian picture. In Section El we apply the method starting 
from an algebraic curve of genus two. 

2. The finite field solution of the discrete KP equation out of 

NONSINGULAR ALGEBRAIC CURVES 

We first briefly recall algebro-geometric construction of solutions of the discrete KP 
equation over finite fields 0011. We discuss in addition a possible degeneracy of the 
linear problem and its consequences. Then we present the Jacobian picture of the 
construction, in which the integrable nature of the equation is evident. We point out 
some aspects of the representation which will help us to construct effectively solutions 
of the equation. 

2.1. General construction. Consider an algebraic projective curve C/K (or simply 
C), absolutely irreducible, nonsingular, of genus g, defined over the finite field K = ¥ q 
with q elements, where q is a power of a prime integer p (see, for example [251 10]). By 
C(K) we denote the set of K-rational points of the curve. By K denote the algebraic 
closure of K, i.e., K = (J^i^Vj an d by C(K) denote the corresponding infinite set 
(often identified with C) of K-rational points of the curve. The action of the Galois 
group G(K/IK) (of automorphisms of K which are identity on K, see []2|) extends 
naturally to action on C(K). 
Let us choose: 

1. four points A e C(K), i = 0,1,2,3, 

2. effective K-rational divisor of order g, i.e., g points fl 7 e C(K), 7 = 1, . . . , g, which 
satisfy the following K-rationality condition 

VaGC7(K/K), a(B y ) = fly. 

As a rule we assume here that all the points used in the construction are distinct and 
in general position. In particular, the divisor ^7=1 -^7 i s non-special. 
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Definition 1. Fix K-rational local parameter t at A . For any integers rii,n 2 ,n 3 G 
Z define the function ifj(ni,n 2 ,n 3 ) as a rational function on the curve C with the 
following properties: 

1. it has pole of the order at most n\ + n 2 + n 3 at Aq, 

2. the first nontrivial coefficient of its expansion in to at Aq is normalized to one, 

3. it has zeros of order at least rii at Ai for i = 1,2, 3, 

4. it has at most simple poles at points £? 7 , j = 1, . . . , g. 

As usual, zero (pole) of a negative order means pole (zero) of the corresponding 
positive order. Correspondingly one should exchange the expressions "at most" and 
"at least" in front of the orders of poles and zeros. By the standard (see e.g., [T]) 
application of the Riemann-Roch theorem (and the general position assumption) 
we conclude that the wave function ip(ni, n 2 , n 3 ) exists and is unique. The function 
ip(ni, n 2 , n 3 ) is K-rational, which follows from K-rationality conditions of sets of points 
in their definition. 

Remark. In what follows we will often normalize functions in a sense of point 2 of 
Definition [TJ 

Fix K-rational local parameters tj at Ai, i = 1,2,3. In the generic case, which we 
assume in the sequel, when the order of the pole of ip{ n i-, n 2i n z) at A is rii + n 2 + n 3 
denote by (l l \ni,n 2 ,n 3 ), i = 0,1,2,3, the K-rational coefficients of expansion of 
il)(ni,n 2 ,n 3 ) at Ai, respectively, i.e., 

^(ni, n 2> n 3 ) = (wi4 i 2+n3) ( l+SfK^,^ ) , 

t \ k=l / 

oo 

V>(ni, n 2 , n 3 ) = tT Cf(m, n 2 , n^, t = 1,2, 3. 

Denote by T; the operator of translation in the variable Ttj, i = 1,2,3, for example 
Tiipijii, n 2 , n 3 ) = ip(ni + l,n 2 ,n 3 ). Uniqueness of the wave function implies the 
following statement. 

Proposition 1. Generically, the function ip satisfies equations 

(1) T^-T,^ + ^H = 0, i^j, i, j = 1,2,3. 

Co 

Remark. When the genericity assumption fails then the linear problem {TJ degener- 
ates, i.e., some of its terms are absent. 

Notice that equation (JTJ gives 

rp Ai) rp >(jf) 

(2) 7§r = -^ ^' = w. 

so So 

Define 

(3) Pl = (-1)^<^^\ 2 = 1,2,3, 
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then equation (j2J) implies existence of a K-valued potential (the r-function) defined 
(up to a multiplicative constant) by formulas 

(4) — = Pi , z = l,2,3. 

T 

Finally, equations (JTJ) give rise to condition 

(5) 1 = 0, 

Pi Pi P2 

which written in terms of the r-function gives the discrete KP equation |T2j called 
also the Hirota equation 

(6) (T\r) (T 2 T 3 r) - (T 2 r) (T 3 T lT ) + (T 3 r) (T\T 2 r) = 0. 

Corollary 2. Equation (JHJ can 6e obtained also from expansion of equation JQ) at 
Afc, where k — 1,2,3, k ^ i, j . 

Remark. Absence of a term in the linear problem (QJ) (see the remark after Proposition 
HJ reflects, due to Corollary in absence of the corresponding term in equation (JHJ). 
This implies that in the non-generic case, when we have not defined the r-function 
yet, we are forced to put it to zero. Let us notice in advance (see the next section) that 
it is in complete analogy with the well known (complex field) interpretation of the 
algebro-geometric solution r of the discrete KP equation as, essentially, the Riemann 
theta function. 

2.2. The Jacobian interpretation. Denote by Div(C) the abelian group of the 
divisors on the curve C and by Pic°(C) the group of eqivalence classes of degree zero 
divisors Div°(C) modulo the principal divisors. There exists [TUEl! an abelian variety 
J{C) of dimension g (the Jacobian of the curve) and an injective map : C —* J{C) 
(the Abel map) such that the extension of to Div(C) establishes an isomorphism 
between Pic°(C) and J(C). Moreover, if there exists a K-rational point A e C(K) of 
the curve, then can be defined by 

C 3 P^0(P) = [P-A]e J(C), 

where [P — A] designates the class of the degree zero divisor P — Am Pic (C). 

Denote by D r (C) the effective divisors of degree r of the curve C and by r the 
extension of to D r (C) 

D r (C) 9Dh r (D) = [D - r ■ A] e J(C). 

The direct image of r is a subvariety W r of dimension r if < r < g, and of 
dimension g if r > g. In particular, W 9 _i defines a divisor in J{C). 

The group Pic°(C; K) of eqivalence classes of K-rational degree zero divisors Div°(C; K) 
modulo the principal K-rational divisors can be identified with the abelian group 
J(C; K) of K-rational points of the Jacobian variety. For finite field K the group 
Pic°(C; K) is finite as well and its order can be found using properties of the zeta 
function of the curve (see, for example |25j). 
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Let us present in this picture the description of the wave function ip and of the r- 
function. We choose the point A as the reference point A and consider the following 
divisor D(ni, n 2 , ^3) G Div°(C;K) of degree zero 

a 

D(nx,n 2 , n 3 ) = m{Ao - M) + n 2 {A - A 2 ) + n 3 {A - A 3 ) + J2 B i ~ 9 ■ A), 

7=1 

with linear dependence on n\, n 2 and n 3 . Its equivalence class in Pic°(C;K) has the 
unique K-rational representant of the form 

a 

X(n h n 2 , n 3 ) = X^ni, n 2 , n 3 ) - g ■ A Q . 

7=1 

This equivalence is given by a function whose divisor reads 

a a 
(7) ni(A a - Aq) + n 2 {A 2 - A ) + n 3 (A 3 - A ) + ^{n u n 2 , n 3 ) - B i ~ °- 

7=1 7=1 

If we normalize such a function at A according to Definition [T] it becomes the 
wave function ip. Notice that some of X 1 could be A which would mean that 
[D(n\, n 2 , n 3 )] G W g _i. This correspondence gives rise to a set of important facts. 

Corollary 3. Evolutions in variables rii, n 2 andn 3 define linear flows in the Jacobian. 

Corollary 4. Points X 1; 7=1, ...,g indicate zeros of the wave function which are 
not specified in the previous construction. 

Corollary 5. // [D(m, n 2 , n 3 )) G W 9 _i then the pole of the wave function at A 
has the order less then (m + n 2 + n 3 ), i.e., we are in the non-generic case, thus 
T(ni,n 2 ,n 3 ) = 0. 

Remark. Notice that because [£>(0,0,0)] G" W g _i then r(0,0,0) ^ 0. 

Remark. If K = C then the algebraic curve C is the compact Riemann surface, the- 
orems of Abel and Jacobi identify the Jacobian with quotient of C 9 by the period 
lattice, and a theorem of Riemann identifies W 9 _i with a certain translate of the zero 
locus of the Riemann theta function (see p]). Then, as we mentioned in remark after 
Corollary O the algebro-geometric solution r of the discrete KP equation becomes, 
with appropriate understanding of its argument via the divisor D(ni, n 2 ,n 3 ) and up 
to a not essential and non-vanishing multiplier, the Riemann theta function (see, 
e.g. [Hi])- I n particular, in such an interpretation the zeros of r are located in points 
of the translate W 9 _i of the Theta divisor. 

Let us discuss periodicity of solutions of the finite field version of the KP equation 
obtained using the above method. Denote by IT, i = 1,2,3, the ranks of cyclic 
subgroups of J(C; K) generated by divisors A; — A , then for arbitrary ki G Z, i = 
1,2,3, 

D{n x + fcilli, n 2 + k 2 U 2 , n 3 + k 3 Il 3 ) ~ £>(m, n 2 , n 3 ). 
In particular, r(ni, n 2 , n 3 ) = implies r{n\ + k\Tli,n 2 + k 2 U 2 ,n 3 + k 3 H 3 ) = 0. 
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There exist unique (normalized at A ) functions hi, i = 1,2,3, with zeros of order 
Ilj at Ai, poles of order IT at A and no other singularities and zeros such that 

(8) ij){n x + feinx, n 2 + k 2 U 2 , n 3 + k 3 U 3 ) = h^hfrhg ^{n x , n 2 , n 3 ). 
Remark. Generalizing above considerations, if for £i G Z, i = 1,2,3, 

3 

J2UA-A )~o, 

i=i 

then (^1,^2,^3) is the period vector of zeros of the r-function and vector of quasi- 
periodicity (in the above sense) of the wave function. 

Equation ((HI) implies quasi-periodicity of the functions Q , i = 1,2,3, 

Co (ni + feilli, n 2 + k 2 ll 2 , n 3 + k 3 ll 3 ) = c^c^c^ Q (ni, n 2 , n 3 ), 
with the (non-zero) factors c^y G K* equal to 

C ( l )j I .SnU, , p=A _ 

The multiplicative group is a cyclic group of order q — 1, therefore the functions 
Co are periodic. Their periods in variable rij are equal to Ilj times the order of the 
subgroup of K* generated by c^j (a divisor of q— 1). Due to possible change of sign 
(see equation Q) the periods of pi can be eventually doubled with respect to the 
corresponding periods of Co ■ Again, periodicity of pi implies quasi-periodicity of r 
with a factor from K*, thus the period of r in variable rii can be maximally q — 1 
times the period of pi in that variable. 

3. A "hyperelliptic" solution of the discrete KP equation 

Our goal here is to demonstrate how does the method described above work. We 
perform all steps of the construction (see also [2] for details) starting from a given 
algebraic curve, which we have chosen to be a hyperelliptic curve, due to relatively 
simple description of Jacobians of such curves [19]. We consider a hyperelliptic curve 
of genus g = 2 but the technical tools used here can be applied directly to hyperelliptic 
curves of arbitrary genus. 

3.1. A hyperelliptic curve and its Jacobian. Consider a hyperelliptic curve C of 
genus g — 2 defined over the field F 7 and given by the equation 

(9) C: v 2 + uv = u 5 + 5u 4 + 6u 2 + u + 3. 

The (u, v) coordinates of its F 7 -rational points are presented in Table [TJ The curve 
has one point at infinity, denoted by 00, whose preimage on the nonsingular model of 
C consists of one point only [21], and where the local uniformizing parameter can be 
chosen as u 2 /v (u is a polynomial function of order 2, and v is a polynomial function 
of order 5). The point opposite (with respect to the hyperelliptic authomorphism) to 
P is denoted by P. The only two special point of the curve are (6, 4) and the infinity 
point 00. 
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P l 







oo 


Po 


1 


(1,1) 


(1,5) 


2 


(2,2) 


(2,3) 


3 


(5,3) 


(5,6) 


4 


(6,4) 


Pa 



Table 1. F 7 -rational points of the curve C. 



i 


Pi 


P l 


Pf 


Pi 


5 


(0,21) 


(0,28) 


P 5 


P 5 


6 


(3,9) 


(3,44) 


Pe 


Pe 


7 


(4, 26) 


(4,33) 


P 7 


Pi 


8 


(7,5) 


(7,44) 


(42,5) 


(42,9) 


9 


(8,22) 


(8,26) 


(43, 29) 


(43, 33) 


10 


(11,5) 


(11,47) 


(46,5) 


(46, 12) 


11 


(12,6) 


(12,45) 


(47,6) 


(47, 10) 


12 


(13,14) 


(13,29) 


(48,35) 


(48, 22) 


13 


(14,8) 


(14,34) 


(35,43) 


(35, 27) 


14 


(15,13) 


(15,28) 


(36,48) 


(36,21) 


15 


(16,17) 


(16,23) 


(37,38) 


(37,30) 


16 


(17,0) 


(17,39) 


(38,0) 


(38, 18) 


17 


(18,4) 


(18,41) 


(39,4) 


(39,20) 


18 


(19,9) 


(19,28) 


(40, 44) 


(40,21) 


19 


(20,12) 


(20,31) 


(41,47) 


(41,24) 


20 


(22,4) 


(22,30) 


(29,4) 


(29,23) 


21 


(25,6) 


(25,32) 


(32,6) 


(32,25) 


22 


(27,7) 


(27,22) 


(34,42) 


(34,29) 



TABLE 2. F 49 -rational points of the curve C (which are not F 7 -rational); 
here P is the opposite of P, and P a denotes its conjugate with respect 
to the lift of the Frobenius automorphism. 



We identify the field F 49 as the extension of F 7 by the polynomial x 2 + 2, i.e., 
F 49 = F 7 [x]/(a; 2 + 2). Let us introduce the following notation: the element k e F 49 
represented by the polynomial j3x + a is denoted by the natural number 7/3 + a. The 
Galois group G(F 4 9/F 7 ) = {id, a}, where a is the Frobenius automorphism, acts on 
elements of F 4 g \ F 7 in the following way: 

F 49 \ F 7 3 k = 7(3 + a i-> a{k) = 7(7 - (3) + a. 

The coordinates of F 49 -rational points of the curve (which are not F 7 -rational) are 
presented in Table [2j 
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In the next step we find the group of the F 7 -rational points J(C; F 7 ) of the Jacobian 
of the curve. The number of its points can be found from the number of F 7 -rational 
and F49-rational points of the curve by application of properties of the zeta function 
of the curve C (see for instance |25, 13;). In our case the curve has 8 F 7 -rational points 
and 74 F4g-rational points which implies the following form of the zeta function ((C; T) 

P(T) 

C(C; T) = (i-r)(i-7T) ' P(T) = 1 + 12T ' + 49T4 ' 

The number # J(C; F 7 ) of the F 7 -rational points points of the Jacobian is equal to 
P(l) = 62, and therefore J(C; F 7 ) is the direct sum of cyclic groups of orders 31 and 
2. 

Let us choose the infinity point oo as the basepoint. The group law in the Jacobian 
of a hyperelliptic curve can be intuitively described in a way which is a higher-genus 
analog of the well known addition operation for points of elliptic curves. We present 
here only its sketch for genus g = 2 and in the generic case of addition of two points 
of J{C) with representants of the form 

Ei = Qi + R { - 2oo, % = 1, 2, 

with all points distinct. If E 3 = Q 3 + R 3 — 2oo is the representant of [E 1 + E 2 ], i.e., 

E 1 + E 2 = (g) + E 3 , 

then 

(10) E x + E 2 + E 3 ~ 0, 

where we have used the fact that for any point P G C(K) of a hyperelliptic curve 
the divisor P + P — 2oo is principal. Therefore, there exists a normalized polynomial 
function / of the order six, thus necessarily of the form 

/ = a + bu + cu 2 + du 3 + ev, 

with divisor given by the left hand side of equation (fTfl|) . Its zeros at Qi and Ri, 
i = 1,2, and the normalization condition allow to fix the coefficients and then to find 
two other zeros Q 3 and R 3 . 

Geometrically, we are looking for two other intersection points of the cubic interpo- 
lating known four points with the hyperelliptic curve. In cases when some points of 
Ei + E 2 are repeated, the interpolation step must be adjusted to ensure tangency to 
the curve with sufficient multiplicity. When divisors have less points then we consider 
the interpolating curve of lower degree (some intersection points are at infinity) . Fi- 
nally, the transition function g is the unique normalized function with the nominator 
equal to / and the denominator being the normalized polynomial function with the 
divisor E 3 + E 3 . 

The full description of the group J(C; F 7 ) is given in Table 03 The divisor D\ = 
Pi — oo generates the subgroup of order 31 and the divisor _D 4 = P 4 — oo generates the 
subgroup of order 2. We present the reduced representants [nD 1 + m_D 4 ] r of elements 
[nDi + mD 4 ] of J(C; F 7 ), where n E {0, 1, ... , 30} and m G {0, 1}. Moreover we write 
down the transition functions g m ( n ) defined by the following divisor equation 

(11) [nDi + mD A ] r + Di = {g m {n)) + [{n + l)D x + mD 4 } r . 
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Also the transition functions for the sums [nDi + mD A ] r + _D 4 can be read off from 
Table El In particular, to find such a transition function (we call it W(0, 1)) with 
n = 29 and m = 0, i.e., 

(12) (1, 5) + (1, 5) + D 4 - 2oo = (W(0, 1)) + (12, 6) + (47, 6) - 2oo, 

we make use of the fact that the analogous transition function for n = 30 and m = 
is 1. Then 

(13) W(0, 1) = ^(29) ■ 1 ■ M29)]- 1 = 2 + 3 " + 4 ^ ± 

o + 4u + M 

where in the last equality we have used equation (JHJ of the curve to get rid of v from 
the denominator. 

3.2. Construction of the wave and r functions. In order to find a solution of 
the discrete KP equation let us fix the following points of the curve C, 

A = oo 1 A 1 = (1,1), A 2 = (2,2), A 3 = (5,3), 

with the uniformizing parameters 

t = u 2 /v, t\ = u — 1, t 2 = u — 2, t 3 = u — 5, 

and 

^ = (12,6), 5 2 = (47,6). 



Because 

A - A 1 ~ 30Di, A - A 2 ~ 21Di + £> 4 , A - A 3 ~ 17Di + D 4 , 
B l + 5 2 - 2A ~ 29^! + L> 4 , 

then the points X\(n\, n%) and X 2 (rii, n 2 , n 3 ) where the wave function ip(ni, n 2 , n 3 ) 
has additional zeros can be found from Table El and from equation 

(14) X 1 (ni,n 2 ,n 3 ) + X 2 (n 1 , n 2 , n 3 ) - 2oo = [nD 1 +mD 4 ] r , 
where n & {0,1, ... , 30} and m G {0, 1} are given by 

(15) n = 29 + 30ni + 21n 2 + 17n 3 ) mod 31, 

(16) m = 1 + n 2 + n 3 mod 2. 

Remark. Notice that the infinity point oo is the Weierstrass point of the curve C, 
which violates the assumption of general position of points used in the construction 
of solutions of the discrete KP equation. This will not change the Jacobian picture of 
the construction but in some situations, which we will point out, will affect uniqueness 
of the wave function. We remark that such a choice is indispensable in reduction of 
the method from the discrete KP equation to the discrete KdV equation (see, for 
example [TK| l3~]). 

To find effectively the wave function we will constraint the range of parameters 
from Z 3 to the parameters of the group of F 7 -rational points of the Jacobian. Let us 
introduce functions h\ and h± corresponding to generators of the two cyclic subgroups 
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n 


[nDi] r 


9o{n) 


[nDi + D4] r 


9i{ n ) 










1 


(6,4) 


— oo 


1 


1 


(1,1) 


— oo 


1 


(1,1) + (6, 4) 


-2oo 


b+bu+3u 2 +v 
6+4u+u 2 


2 


(1,1) + (1,1) 


- 2oo 


u+5u 2 +v 


(12,45) + (47,10) 


-2oo 


l+5u 2 +u 


(2+uj 2 


2+5u+u 2 


3 


(5, 6) + (5, 6) 


- 2oo 


l+u+4u 2 +u 
(2+u) (5+u) 


(15,28) + (36,21) 


-2oo 


6u 2 +v 
2+u 2 


4 


(2, 3) + (5, 3) 


- 2oo 


2+4u 2 +t> 


(7, 44) + (42,9) 


- 2oo 


5+U+D 


5+4U+U 2 


4+6u+u 2 


5 


(19,9) + (40,44) 


- 2oo 


4u+2u 2 +t> 


(11,5) + (46,5) 


- 2oo 


6+6u+u 2 +u 


5+5u+u 2 


3+6u+u 2 


6 


(22, 4) + (29, 4) 


- 2oo 


5+2u+6u 2 +u 


(18,41) + (39,20) 


-2oo 


5+3u+5u 2 +v 


(2+u) (5+u) 


5+3u+u 2 


7 


(2, 3) + (5, 6) 


- 2oo 


5+6u+2u 2 +u 
5+2u+u 2 


(16,17) + (37,38) 


-2oo 


5+4u+4u 2 +u 
3+u+u 2 


8 


(27,22) + (34, 29) 


- 2oo 


l+3u+2u 2 +u 
1+u 2 


(17,39) + (38,18) 


-2oo 


3+2u+u 2 +u 
(5+u)(6+u) 


9 


(14,34) + (35,27) 


- 2oo 


l + 5u+l> 


(1,5) + (2, 2) 


- 2oo 


6 + u 


(l+u)(5+u) 


10 


(2, 2) + (6, 4) 


- 2oo 


3+5u+5u 2 +u 
(5+u) 2 


(2,2). 


— oo 


1 


11 


(2, 3) + (2, 3) 


- 2oo 


6+u+6u 2 +v 


(1,1) + (2, 2) 


- 2oo 


4+2u 2 +u 


3+2u+u 2 


3+5u+u 2 


12 


(13, 14) + (48, 35) 


- 2oo 


3+6u+4u 2 +u 


(8,22) + (43,29) 


-2oo 


2+4u+u 


2+2u+u 2 


(2+u)(6+u) 


13 


(20,12) + (41,47) 


- 2oo 


5u+u 2 +v 


(1,5) + (5, 3) 


- 2oo 


6 + u 


(l+u)(2+u) 


14 


(5, 3) + (6, 4) 


- 2oo 


6+5u+2u 2 +u 


(5,3)- 


— oo 


1 


6+6u+u 2 


15 


(25,32) + (32,25) 


- 2oo 


5+u 2 +v 


(1,1) + (5, 3) 


- 2oo 


u+5u 2 +v 


6+6U+U 2 


(2+u)(6+u) 


16 


(25,6) + (32,6) 


- 2oo 


6+5u+2u 2 +u 
(1+u) (2+u) 


(1,5) + (5, 6) 


- 2oo 


+ U 


17 


(5, 6) + (6, 4) 


- 2oo 


9 

5u+u +u 

2+2U+U 2 


(5,6) 


— oo 


1 


18 


(20,31) + (41,24) 


- 2oo 


3+6u+4u 2 +u 
3+2u+u 2 


(1,1) + (5, 6) 


- 2oo 


2+4u+u 
3+5u+u 2 


19 


(13,29) + (48,22) 


- 2oo 


6+u+6u 2 +v 


(8,26) + (43,33) 


-2oo 


4+2u 2 +u 


(5+u) 2 


(5+u)(6+u) 


20 


(2,2) + (2,2) 


- 2oo 


3+5u+5u 2 +v 

(1+u) (5+u) 
l+5u+t) 

1+u 2 


(1,5) + (2, 3) 


- 2oo 


6 + u 


21 


(2, 3) + (6, 4) 


- 2oo 


(2,3) 


— oo 


1 

± 


22 


(14,8) + (35,43) 


- 2oo 


l+3u+2u 2 +u 


(1,1) + (2, 3) 


- 2oo 


3+2u+u 2 +u 


5+2u+u 2 


3+u+u 2 


23 


(27, 7) + (34,42) 


- 2oo 


5+6u+2u 2 +u 


(17,0) + (38,0) 


- 2oo 


5+4u+4u 2 +D 


(2+u)(5+u) 


5+3u+u 2 


24 


(2, 2) + (5, 3) 


- 2oo 


5+2u+6u 2 +u 
5+5u+u 2 
4u+2u 2 +i> 


(16,23) + (37,30) 


-2oo 


b+3u+5u 2 +v 
3+6u+u 2 


25 


(22,30) + (29,23) 


- 2oo 


(18,4) + (39,4) 


- 2oo 


6+6u+u 2 +u 


5+4u+u 2 
2+4u 2 +v 
(2+u)(5+u) 


4+6u+u 2 


26 


(19, 28) + (40, 21) 


- 2oo 


(11,47) + (46, 12) 


-2oo 


5+u+v 
2+u 2 


27 


(2, 2) + (5, 6) 


- 2oo 


l+u+4u 2 +u 


(7, 5) + (42, 5) 


- 2oo 


6u 2 +v 


(2+u) 2 


2+5u+u 2 


28 


(5, 3) + (5, 3) 


- 2oo 


u+5u 2 +v 


(15,13) + (36,48) 


-2oo 


l+5u 2 +u 


(6+u) 2 


6+4u+u 2 


29 


(1,5) + (1,5) 


- 2oo 


(6 + u) 


(12, 6) + (47,6) • 


- 2oo 


5+5u+3u 2 +v 
(l+u)(6+u) 


30 


(1,5) 


— oo 


(6 + u) 


(1,5) + (6, 4) 


- 2oo 


6 + U 



Table 3. The group J(C;F 7 ) of F 7 -rational points of the Jacobian as 
the simple sum of its cyclic subgroups. 
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of J(C; F 7 ). The function hi with the divisor 31-Di and normalized at the infinity point 
is equal to the product Yli=o9o^) an( ^ rea ds 

hi = 1 + 2u + u 2 + Au 3 + 3u 5 + m 6 + 3m 7 + u 8 + 4w 9 + 

4u 10 + 2m 11 + 5m 12 + 2u n + 4u 14 + 3m 15 + [hu + 2u 2 + 

5w 3 + 4m 5 + 6m 6 + 4m 7 + 3m 9 + 5m 10 + 5m 11 + 4m 12 + u 13 )v, 

where we also used equation of the curve Q to reduce higher order terms in v. The 
normalized function h^ with the divisor 2D 4 is 

hi = u — 6. 

Let us introduce other auxiliary functions f 2 and f 3 to factorize the zeros of the wave 
function at A 2 and A 3 . Notice that 

A 2 + 21D 1 + £> 4 - 00 ~ 0, 

which implies that there exists a polynomial function on C with simple zero at A 2 and 
other zeros in the distinguished (by our choice of description of J(C; F 7 )) points (1, 1) 
and (6,4). Define f 2 as the unique such function normalized at the infinity point 00, 
then 

f 2 = 1 + 5u + u 2 + Au A + 6u 5 + 4u 6 + Au 7 + 3u 8 + Au 9 + 

Qu 11 + (6 + 4m + 2m 2 + 5m 3 + 6m 4 + Qu 6 + u 7 + u 8 + u 9 )v. 
Similarly we define the normalized function 

/ 3 = 1 + Q U + 2m 2 + 6m 5 + u 6 + 5m 7 + 5m 8 + 4m 9 + (4 + 3m + 5m 2 + 4m 5 + 2m 6 + u 7 )v, 

with the divisor A 3 + 17 Di + _D 4 — 00. 

Uniqueness of the wave function ip implies that it can be decomposed as follows 

(17) ip(n 1} n 2 ,n 3 ) = W(mi,m 2 ), 

where W / (m 1 ,m 2 ) is the unique normalized function with the divisor 

(18) miDi + m 2 D 4 + Yi(m u m 2 ) + Y 2 (m u m 2 ) - (12, 6) - (47, 6), 
where 

(19) Yi(mi,m 2 ) + Y 2 (m 1 , m 2 ) = Xi(m, n 2 , n 3 ) + X 2 (rii, n 2 , n 3 ), 
and the new variables mi i m 2 are given by 

(20) 21t7 2 + 17n 3 - ni = 31p- mi, m x G {0, 1, . . . , 30}, 

(21) n 2 + n 3 = 2q-m 2l m 2 G {0, 1}. 

To find the functions W(mi, m 2 ) for all mi 6 {0,1,..., 30} and m 2 G {0, 1} let us 
first notice that W(0, 0) = 1 and W(0, 1) is indeed the function found in equation 
(fl3|) . For mi G {1, ... , 30} and m 2 G {0, 1} define the functions w m2 {mi) as follows 

W(mi,m 2 ) = w rri2 (m 1 )W(mi - l,m 2 ). 
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Equations JUJ, (JHJ-JTBJ) and (fT%l) - (|2*Tl) imply that for such range of mi and m 2 we 

have 

g m (n) = w m2 (mi), 

where 

m 2 = 1 — m mod 2, mi = 29 — n mod 31. 
Finally, under identification w m2 (0) = W(0,m 2 ) we obtain 

mi 

W(mi,m 2 ) = [jw^fi), 

which, together with factorization fTTjl . gives the wave function if) for all values of 

(ni,n 2 ,n 3 ) G Z 3 . 

Remark. For (mi,m 2 ) = (29, 1) we have Xi = X 2 = oo. Because the infinity point 
oo is the Weierstrass point of order two, there exist functions with divisor of poles 
equal to 2oo. This means that if) is not uniqely determined in this case. However it 
is natural to keep the divisor of if>, and therefore if) itself, exactly like it is given from 
the flow on Jacobian. Notice that because for X\ = X 2 = oo we stay in the divisor 
W 9 _i then this ambiguity does not affect construction of the r-function. 

The coefficients Co (^1,^2,^3)1 k = 1,2,3, of expansion of the wave function can 
be obtained from factorization (JTYj) and are given by 

(22) tfV,n 2 ,n 3 ) = 6^4* (^|^) 

(23) Ci 2) (ni,n 2 ,n 3 ) = 6^+H p W{m Xl m 2 ) \ t2=0 , 

(24) Ci 3) (ni,n 2 ,n 3 ) = 5 nM 3 p 6 q W(m u m 2 ) | t3=0 . 

Using definition of the r function for nonzero p^, i.e. equation (^J), and putting r = 
for points of the divisor W 9 _i we obtain the corresponding F 7 - valued solution of the 
discrete KP equation (0). This r-function is presented in Figure [TJ Notice that due 
to quasi-periodicity of the r-function, we have to calculate the solution in this way 
only for a finite range of values of the variables. 

The periods IT, i = 1,2,3, of zeros of the r-function are, respectively, 31, 62 and 
62. Equivalently, the "period vectors" of zeros can be choosen as 

Because C(i)i = 6, C( 2 )i = 3 and C( 3 )i = 5 then periods of the functions pi, i = 1,2,3, 
in that variable are, respectively, 2 • 31, 3-31-2 and 6 • 31. Moreover, we have 
r(31, 0, 0) = 3 which gives r(62, 0, 0) = 5 = 3 2 • 6 31 mod 7, and therefore the period 
of r in rii is 6 • 2 • 31. 
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□ -0, B-l, -2, -3, B-4, -5, -6. 

Figure 1. F 7 -valued solutions of discrete KP equation out of genus 
two hiperelliptic curve C. Variables n x (directed to the right) and n 2 
(directed up) take values from to 34. Subsequent figures are for values 
of n 3 = -1,0, 1,2. 
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